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摘要 :利用 大 气 和 水 文 模型 定量 描述 陆 表 相关 变量 变化 规律 一 直 是 大 气 科 学 和 水 文学 界 的 研究 热点 。 狼 而 ,由 于 我 国 西部 地 区 
站 点 匮乏 ,传统 气象 观测 站 点 已 不 能 满足 大 尺度 地 表 分 量 高 精度 模拟 分 析 的 需求 。 建 立 SWAT 模型 中 国 大 气 同 化 驱动 数据 集 
(China Meteorological Assimilation Driving Datasets for the SWAT model, CMADS) 驱动 SWAT(Soiláànd Water Assessment Tool) fi 
型 (简称 为 CMADS+SWAT 模式 ) ,选取 传统 气象 站 点 稀缺 的 新 疆 精 博 河流 域 为 并 区 , 兄 成 流域 各 地 表 分 量 (如 土壤 湿度 . 雪 深 、 
融雪 ) 校准 .验证 及 其 时 空 关 系 提取 与 分 析 。 分 析 发 现 :CMADS 数据 集 可 很 好 地 驱动 KE SWAT 模式 完成 本 地 化 工作 。 其 
中 ,CMADS+SWAT 模式 在 月 尺度 上 总 体 NSE 效率 系数 均 在 0.659—0.942, 目 尺 度 也 均 在 0.526 一 0.815。 对 流域 内 土壤 湿度 和 
融雪 过 程 进 行 相关 分 析 发 现 : 精 博 河流 域 土壤 湿度 在 年 内 3 一 4 月 份 达到 其 第 宇 次 峰值 ,主要 贡献 来 自 于 流域 内 高 山 融 雪 现 象 ; 
融雪 期 结束 后 ,流域 降水 量 增 加 ,伴随 气温 上 升 等 现象 导致 土壤 温度 旺 现 波动 态势 ,至 10 月 中 旬 冷 空气 过 境 产 生 较 大 降水 
( 雪 ) ,最 终 使 土壤 水 转变 为 冻 土 ,直至 次 年 接近 融雪 期 ,土壤 水 再 次 增加 且 到 融雪 过 程 结束 。 一 方面 证 明 CMADS+SWAT 模式 
可 有 效 提高 SWAT 水 文 模型 在 我 国 西北 干旱 区 (站 点 稀缺 区 域 ) 的 表现 能 力 ,为 一 方面 理 清 了 精 博 河流 域 相关 地 表 分 量 ( 土壤 
湿度 .蒸发 等 ) 时 空 演 变 规律 。 本 研究 对 我 国 大 气 水 文学 种 发 展 将 起 到 一 定 的 科学 促进 作用 。 
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Abstract: Describing the changing rules of land surface variables by using meteorological and hydrological models has 
always been/@ research hotspot in the fields of atmospheric and hydrological science. However, due to the scarcity of weather 
stations, in West China, traditional weather stations cannot satisfy the requirements of high precision simulation of large scale 
land surface, components. Using the Jing and Bo River Basin in Xinjiang as a research area, the present study used the 
China Meteorological Assimilation Driving Datasets for the SWAT model ( CMADS) to drive the Soil and Water Assessment 
Tool. (SWAT) model, and then completed calibration, verification, and time-space relation analysis of each land surface 
component (such as soil moisture content, snow depth, and snow melt). Our analysis showed that the CMADS dataset can 


drive and calibrate the SWAT mode land complete localization work in Jin and Bo River Basin. The NSE efficiency 
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coefficients of the SWAT model driven by the CMADS dataset were generally controlled between 0.659 and 0.942 at the 
monthly scale, and were also controlled between 0.526 and 0.815 at the daily scale. Furthermore, our analysis indicated 
that the soil moisture content would reach a high level for the first time between March and April each year, which 1s mainly 
caused by snow melting in the high mountains. However, when snow melting finished, due to an increase in precipitation 
and temperature, soil temperature fluctuated until the middle of October, when cold air brought considerable precipitation 
and snow. Finally, soil water was transformed into frozen soil until the snow melting period in the following year. Thereafter, 
soil water would increase again until the end of the snow melting period. On one hand, this study verified that the CMA DS+ 
SWAT mode can enhance the performance ability of the SWAT model in the arid areas of northwestern China, which lacks 
weather stations. On the other hand, the study provided a scientific explanation for time-space changing rules of land surface 
components ( such as soil moisture and evaporation) in the Jing and Bo River Basin. The findings of this research will play a 


certain role in promoting the development of China's meteorological and hydrological sciences. 


Key Words: CMADS; SWAT; Jing and Bo River Basin; surface component 


近年 来 ,全球 范围 内 均 面 临 严 重 的 生态 与 水 危机 ,干旱 区 陆 面 过 程 变化 无 论 对 整个 内 陆 河 流水 循环 亦 或 
是 对 整个 干旱 区 生态 环境 演变 均 具 重要 影响 。 其 中 ,干旱 内 陆 区 水 量 这 k 山 加 速 又 减 与 恶化 也 促使 该 地 区 突 
发 性 水 事件 频繁 发 生 。 干旱 区 绿洲 构造 独特 ,水 与 能 量 循环 具备 其 自身 的 特殊 规律 ,系统 分 析 新 疆 干 旱 区 地 
表 水 资源 相关 分 量 时 空 变化 意义 重大 ,这 将 为 我 国 实现 内 陆 干 旱 区 流域 生态 水 文 可 持续 发 展 提供 重要 科学 依 
据 。 精 博 河 流域 位 于 新 疆 地 区 西部 ,封闭 流域 总 面积 达 25 万 km*。 流 域 在 西北, 南 等 3 个 方向 均 为 高 山 ,中 
间 为 谷地 ,东边 为 例 地 。 地 貌 及 地 形 宏观 上 与 准噶尔 倪 地 连 为 = 条 ,大 体 分 为 山谷 .分 三 大 地 貌 。 流 域 交汇 
处 正好 位 于 阿拉 山口 大 风 下 沉 通 道 , 故 大 风灾 害 是 该 流 虞 最 大 的 气象 灾害 。 此 外 ,由 于 该 流域 生态 结构 极其 
脆弱 ,导致 研究 区 内 风沙 频繁 . 盐 渍 化 及 攻 漠 化 现象 严重 。 近 年 来 ,由 于 气候 变 暖 加 剧 及 人 类 垦荒 活动 不 断 ， 
导致 流域 内 艾 比 湖人 湖水 量 不 断 减少 ,生态 环境 持续 屎 化 。 据 不 完全 统计 ,日 20 世纪 50 ERE 70 年 代 以 
来 , 精 博 河 流域 内 芯 比 湖面 积 萎 缩 了 近 678km“ 哆 比 湖 储 水 量 共 损耗 近 23 亿 m 。 水 资源 的 急剧 下 降 导 致 艾 
比 湖水 体 矿 化 度 急剧 攀升 ,湿地 面积 嗓 然 减 小 ,湖泊 对 气候 的 调控 作用 降低 ,对 流域 内 生态 系统 及 当地 居民 的 
生存 环境 造成 了 极 大 的 破坏 影响 4 并 最 终 导 致 精 博 河 生态 环境 迅速 恶化 ,发 生 了 生物 链 断 层 多样 性 减少 ,元 
省 化 过 程 加 剧 等 现象 。 由 于 精 博 河流 域 对 新 疆 生 态 平 衡 及 社会 经 济 发 展 均 具 有 重大 页 献 ,利用 高 精度 大 气 数 
据 驱 动 成 熟 水 文 模型 的 输 乔 结果 可 为 流域 提供 一 定 的 可 持续 发 展 决 案 , 促 进 精 博 河 流域 地 表 过 程 分 量 ( 和 后 
态 、 水 文 等 ) 合理 发 展 。 

由 于 精 博 河流 域 传 统 和 气象 站 点 稀缺 ,流域 下 热 面 分 异 较 大 ,加 之 气候 变化 及 人 为 活动 的 影响 ,增强 了 该 区 
域 模 式 输出 与 评估 的 不 确定 性 。 众 多 研究 证 明 : 融 合 更 多 观测 数据 的 大 气 同 化 数据 集 可 更 大 程度 提高 模 
式 输出 精度 坚 3 。 当 前 ,国内 外 大 气 再 分 析 种 类 繁多 ,如 : NCEP 的 NCAR/DOE'*". ,欧洲 预报 中 心 ERA 一 15、 
ERA—40.,ERA--Interim 再 分 析 资 料 551 日 本 JRA—25 再 分 析 资 料 ' 引 及 陆 面 模式 驱动 场 Princeton 数据 集 。 由 
于 传统 火气 观测 人 站 并 未 布 满 全 球 , 以 上 再 分 析 数 据 集 对 气象 缺 测 区 气候 场景 分 析 及 应 用 提供 了 重要 数据 支 
BE. 然而 ,对 于 区 域 精细 化 模拟 而 言 ,以 上 再 分 析 数 据 的 准确 性 示 待 提高 。Pal 等 … 利用 区 域 气候 模式 
RegCM3 对 东亚 季风 区 冬 . 夏 两 个 季节 内 降水 量 的 月 变化 规律 在 季节 尺度 上 进行 了 模拟 与 评估 后 发 现 ; 区 域 
模式 RegCM3 对 降水 量 存在 较 大 误差 ,这 种 现象 在 冬季 显得 尤为 凸显 。 赵 天 保 等 利用 NCEP ERA 及 JRA 等 
数据 集 的 地 面 压强 及 温度 要 素 在 中 国 区 域 可 用 性 进行 评估 后 发 现 :3 种 再 分 析 数 据 集 存 在 明显 的 季节 性 和 区 
域 性 差异 中。 施 晓 上 晕 等 “利用 各 类 技术 手段 将 NCAR 驱动 数据 (包括 :空气 温度 风速) 在 中 国 区 域 的 适 
用 性 进行 评估 后 发 现 : 风 场 数据 距 平 误差 与 海拔 表现 为 负 相 关 现 象 。 以 上 人 研究 均 表 明 . 各 类 再 分 析 产 品 虽 可 
较 好 地 反映 出 大 尺度 气象 要 素 空 间 模 态 ,然而 ,由 于 大 部 分 再 分 析 数 据 集 并 未 利用 我 国 全 部 区 域 自动 站 进行 
订正 , 故 不 能 较 好 地 反映 我 国 近 地 表 真实 气象 状况 。 综 上 所 述 : 利 用 更 接近 我 国 真实 气象 场 的 大 气 数据 集 驱 


http ://www.ecologica.cn 


201711.02248v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


7116 生 态 学 报 37 卷 


动 成 熟 的 水 文 模型 ,可 更 好 的 对 我 国 大 尺度 地 表 分 量 时 空 演变 规律 进行 细致 化 分 析 | 。 

本 人 研究 以 站 点 稀缺 .下 垫 面 分 异 较 大 的 新 疆 博 精 河流 域 作为 典型 验证 及 分 析 区 ,利用 SWAT 模型 中 国 大 
气 同化 驱动 数据 集 (CMADS) ”驱动 在 我 国 应 用 普遍 的 SW AT 水 文 模型 ,以 期 在 为 该 流域 提供 有 效 生态 水 文 
系统 恢复 策略 的 同时 ,验证 CMADS 数据 集 在 中 国 西部 地 区 的 适用 性 。 


1 研究 区 介绍 


精 博 河流 域 地 处 西天 山北 坡 ( 图 1) , 4 F 81946'—83?51'E,44?02'—45?10'N 之 间 。 总 面积 达 11275km?, 
流域 内 降水 多 为 来 目 北 冰 洋 和 大 西洋 的 水 汽 ,呈现 总 体 山区 较 平 原 多 ,西部 较 东 部 多 , 阴 坡 较 阳 坡 多 的 态势 。 
精 河 和 博 尔 塔 拉 地 区 的 冰川 面积 均 分 别 为 96.2km* 和 110.3km*。 其 中 , 精 河 冰川 补给 量 在 年 际 考 约 为 0.96 亿 
m , 博 尔 塔 拉 河 约 为 1.05 亿 m ,分 别 占据 了 其 本 身 河川 总 径流 的 20.6% 和 21.4%。 近 年 来 ,流域 内 点 比 湖 湖 
水 矿 化 度 在 逐年 升 高 ,流域 生态 与 生活 用 水 矛盾 也 随 之 日 益 凸 显 。 流 域 土地 退化 现象 严重 ,其 中 ;区 比 湖 流 域 
湖底 已 有 1500km ”区域 退化 沦 为 盐 漠 , 盐 渍 化 面积 达 71km 。 
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1 博 、 精 河流 域 示 意图 


Fig.1 Schematic diagram of the Jing and Bo River Basin 


Yr A3 , "r8 DORT RS E387 i398 77 E 8 ^p 39 88 DC BAL EA UE 2 DU T 8 2S BEA EOS (3X DORIA, ) 及 
ERARI. HTAA 4 个 气象 站 ,导致 对 该 区 域 还 未 进行 较为 合理 可 信 的 水 文 过 程 模拟 与 分 
析 工 作 。 本 研究 将 基于 面 回 东亚 的 高 分 辨 卒 大 气 驱 动 场 CMADS 数据 集 驱 动 SWAT 模型 ,对 精 博 河水 文 相关 
分 量 进行 细致 模拟 ,以 期 系统 分 析 精 博 河流 域 各 地 表 要 素 时 空 演变 特征 ,同时 为 我 国 高 分 辩 率 气象 数据 驱动 
复杂 救 广 模 式 研 究 打 下 坚实 基础 。 


2. 材料 与 方法 


2.1 SWAT 模型 地 表 输 入 数据 
2.1.1 数字 高 程 模型 数据 

SWAT 模型 输入 数据 主要 包括 数字 高 程 模型 (DEM ) 土壤 数据 、 土 地 利用 数据 及 气象 驱动 场 。 其 中 , 数 
字 高 程 模 型 选取 自 SRTM —(90m) DEM( FÉ 2) , 该 数据 来 源 于 CGIAR—CSI SRTM 90m 高 程 数 据 库 ( http:// 
srtm.csi.cgiar.org/ SELECTION/ inputCoord.asp ) 。 
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2 研究 区 高 程 及 SWAT 模型 其 他 输入 数据 
Fig.2 The elevation of the study area and the input data for SWAT model 


2.1.2 ” 精 博 河流 域 土壤 分 布 

土壤 物理 特性 决定 了 SWAT 模型 不 同 水 文 响应 单元 的 产 时 世 流 特征 ,同时 也 是 SWAT 模型 定义 水 文 啊 应 
单元 (HRU ) 的 参考 依据 。 本 研究 选用 的 土壤 输入 数据 内 基于 世界 土壤 数据 库 (HWSD ) 的 中 国土 壤 数 据 集 
(V1.1) (图 3A)'”。 土 地 利用 来 源 于 精 博 河流 域 管理 处 ,土地 利用 数据 又 加 了 我 国 第 二 次 冰川 数据 集 "” Lee 
加 后 的 土地 利用 分 布 情况 可 参见 图 3B。 为 保证 模式 分 辨 率 一 致 ,本 人 研究 将 DEM .土壤 与 土地 利用 等 数据 空 
间 分 状 率 均 统一 为 Tm ,平面 投影 坐标 统 生 设 顽 为 WGS_1984_UTM _Zone_44N。 
2.1.3 ”大气 驱动 数据 

本 研究 使 用 的 大 气 输入 数据 SWAT 模型 中 国 大 气 同 化 驱动 数据 集 ( China Meteorological Assimilation 
Driving Datasets for the SWAT imodel, CMADS) ,该 数据 集 引 入 了 中 国 陆 面 同化 系统 ,保证 了 气象 输入 数据 的 可 
HEE ”。CMADS 站 点 具体 信息 克 表 1 KK 2, 


x1 CMADS 数据 集 信 息 
Table 1 The information of CMADS driving data 


数据 集 CMADS ( V1.0) 
Datasets China Meteorological Assimilation Diving Datasets for the SWAT model ( Version 1.0) 
提供 要 素 Provide elements 日 最 高 / 低 气 温 .日 平均 相对 湿度 .日 均 风 速 .日 累计 降水 量 .日 累计 太阳 辐射 
原始 空间 范围 Original space range 0—65?0'N,60—160?0'E 
本 研究 空间 范围 /Spatialgcope 81°46'—83°51'E  44202'—45?10'N 
Ehi NL E Time scale 2008.1.1—2014.12.31 
原始 分 辨 率 Original resolution 1/3?,1/4?,1/8?,1/16? 
本 人 研究 使 用 空间 分 辨 率 Spatial resolution in this study 1/3? 
站 点 数 Station number 23 个 
CMADS V1.0 版 本 空间 分 辨 率 为 1/3 度 , 时 间 分 辨 率 为 逐日 ,数据 尺度 为 2008 一 2014 年 。SWAT 模型 在 


精 博 河流 域 共 读 取 了 CMADS 数据 集 近 23 个 站 点 要 素 ( 温 湿 、 风 、 降 水 及 辐射 数据 ) 。 其 中 ,CMADS 气象 站 
点 具体 信息 (站 点 名 称 及 经 纬度 ) 如 表 2 所 示 。 
2.1.4 水 文 验 证 数据 

本 研究 将 精 博 河流 域内 水 文 站 逐日 流量 资料 作为 模型 率 定 及 验证 资料 , 台 站 信息 如 表 3 所 示 。 
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图 3 精 博 河流 域 土壤 分 布 及 精 博 河流 域 土地 利用 分 布 
Fig.3 Soil distribution in the Jing and Bo River Basin, and Land use distribution in the Jing and Bo River Basin 
iE :SWGR : 细 芽 草 , Slender Wheatgrass ; PAST ; 1437; , Pasture; HAY ; TF , Hay; WETL : HA IEH , Wetlands-Mixed ; WATR : 7K 3 , Water; AGRL ; H 
他 作物 , Agricultural Land-Generic ; AGRR ; "P SEJEJ , Agricultural Land-Row Crops; URML :居民 区 , Residential-Med/Low Density; UIDU ; 工业 用 
Ji , Industrial ; UTRN :交通 用 地 ,Transportation; FRST ; 林地 , Forest-Mixed; RNGE ; 草场 , Range-Grasses; BARO ; ZR Z1, Bare rock; DESE: 荒漠 ， 
Desert; SALA :盐碱地 ,Saline land; GLSN :冰川 /永久 积 雪 ,Glacier and snow; FRSD :落叶 林 , Forest-Deciduous ; PINE ; £| HER, Pine; RNGB :灌木 
Jk, Range-Brush ; ORCD :果园 , Orchard ; FRSE : 常 绿林 ,Forest-Evergreen 
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表 2 精 博 河流 域 CMADS 站 点 分 布 
Table 2 The station distribution of CMADS in Jing and Bo River Basin 


eet 纬度 dm 经 度 EC 高 程 /m eke 纬度 NC) 经 度 高 程 /m 
of CMADS Latitude Longitude Elevation of CMADS Latitude Longitude Elevation 
136—62 44.99 80.34 2968.00 136—68 44.99 82.34 583.00 
136—65 44.99 81.34 1074.00 135—67 44.65 82.01 504.00 
136—69 44.99 82.68 184.00 135—68 44.65 82.34 296.00 
137—967 45.32 82.01 1863.00 135—64 44.65 81.01 2072.00 
136—66 44.99 81.68 735.00 135—606 44.65 81.68 1432.00 
136—71 44.99 83.34 232.00 134—69 44.32 82.68 936.00 
136—967 44.99 82.01 649.00 134—68 44.32 82.34 1652.00 
136—70 44.99 83.01 206.00 134—67 44.32 82.01 2081.00 
136—72 44.99 83.67 249.00 134—70 44.32 83.01 1898.00 
135—71 44.65 83.34 236.00 134—71 44.32 83.34 3057.00 
135—70 44.65 83.01 206.00 133—70 43.99 83.01 2338.00 

135—060 44.65 82.68 263.00 


CMADS 站 点 名 为 CMADSV1.0 数据 集 在 其 空间 范围 (东亚 区 域 ) 内 唯一 标识 符 


表 3 精 博 河流 域 水 文 站 相关 信息 


Table 3 Hydrological stations data statistics in Jing and Bo River Basin 


水 文 站 纬度 经 度 台 站 海拔 /]m 数据 时 段 /a 
Hydrological station Latitude Longitude Station elevation Data interval 
ili 5& Wenquan 44*59'N 81°02% 1310 2009—2013 
精 河 Jinghe 44?22'N 82?*55'E 620 2009—2013 


2.0 SWAT 水 文 模式 介绍 

SWAT 模型 概 化 的 将 相同 土地 柳 被 利 儿 s 二 二 分 布 及 管理 方式 的 均一 性 单元 视 为 同一 类 水 文 响应 单元 
(HRU) ,模型 认为 每 一 种 HRU 内 所 有 水量 平衡 过 程 具 一 致 性 特征 。SWAT 模型 提供 了 3 种 潜在 蒸发 计算 
方法 , 主要 包括 Penman-Monteith-34; 7", Priestley-Taylor 法 及 Hargreaves 法 ~ 。 在 本 研究 中 , 由 于 
CMADS 数据 集 可 提供 全 部 输入 要 素 (而 传统 气象 站 由 于 无 法 提供 太阳 种 射 数据 仅 能 利用 天 气 发 生 器 模拟 )， 
本 研究 选取 Penman-Monteith 活 帮 为 蒸发 模拟 的 方法 ,以 最 大 可 能 提高 模式 异 拟 精度 52 。 


3 SWAT 模型 方案 设置 


本 研究 区 流域 面积 经 过 ArcGIS 提取 后 为 2482366.4052hm2 , 共 划 分 39 个 子 流 域 及 1648 个 水 文 响应 单 
元 。 由 于 降水 输入 数据 为 逐日 数据 ,地 表 径流 模拟 方法 选取 径流 曲线 数 模型 计算 ,地 表 径 流 将 被 分 别 在 不 同 
水 文 啊 雇 单元 电 进 行 模拟 ,最 终 汇流 到 主 河道 。 最 后 ,选取 基于 连续 方程 的 河 槽 苹 量 法 计算 主 河道 水 量 演进 。 
SW AT d ARIETE Centriod 法 将 流域 内 单 点 气象 数据 空间 离散 插值 到 整个 流域 ,为 减少 空间 离散 插值 之 来 的 误 
25 (JU Be SEA DC) ,本 人 研究 针对 精 博 河流 域 划分 出 多 个 高 程 沉 以 期 精确 识别 不 同 高 程 囊 降水 量 分 布 状 态 。 

综合 CMADS 驱动 场 时 间 尺 度 (2008 一 2014 年 ) 及 径流 观测 数据 的 重合 区 间 ,同时 为 使 模拟 初期 所 有 水 文 
过 程 从 初始 状态 进入 平衡 状态 ,本 研究 将 预 热 期 设 定 为 1 年 ( 即 2008 年 ) ,校准 期 设 为 2009 一 2010 年 ,验证 期 
设 定 为 2011 一 2013 年 。 
3.1 敏感 性 分 析 

本 人 研究 利用 SWAT-CUP 对 SWAT 模式 进行 模式 率 定 工作 。SWAT-CUP 是 EWAGE 研究 所 针对 SWAT Bi 
型 研发 的 一 款 自 动 校准 与 不 确定 性 分 析 的 程序 。 在 本 研究 中 ,重点 对 12 个 与 径流 相关 的 参数 进行 参数 率 
E ,并 得 到 模式 参数 最 终 取 值 ( 表 4) 。 
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3.2 ”模型 校准 

在 模型 率 定期 (2009 一 2010 年 ) 对 模型 参数 进行 校准 工作 。 笔 者 先 利 用 SWAT-CUP 率 定 软 件 对 精 博 河 
流域 两 个 测 站 同时 进行 参数 校准 ,校准 过 程 中 兼顾 参考 流域 年 均 莱 发 量 及 降水 量 情况 下 进行 宏观 调整 ,首先 
考虑 年 雁 发 量 和 径流 量 的 关系 比 ,保证 模拟 结 末 符合 精 博 河流 域 年 总 燕 发 量 、 降 水 量 及 径流 量 接近 实际 真实 
情况 。 分 然后 延展 至 月 ,在 月 数据 校准 准确 后 进行 日 参数 率 定 及 微调 。 分 析 发 现 :SWAT 模型 降水 递减 率 
(PLAPS) 参 数 最 优 取 值 为 44.5mm/km。 温 度 梯度 (TLAPS) 最 优 值 为 -4.3039%C/km( 表 4) ,与 流域 内 多 年 实 
味 均 值 吻合 。 


表 4 SWAT 模型 参数 最 终 取 值 结果 
Table 4 The final results of the model parameters within SWAT model 


变量 名 参数 定义 最 终 优 选 值 
Variable Parameter definition Final optimal value 
CN2.mgt SCS 径流 曲线 值 34 
ALPHA, BF.gw 基 流 a BIF (days) 0.453 
GW_DELAY.gw 地 下 水 补给 延迟 时 间 (days ) 42 
GWQMN.gw 地 下 水 汇 和 人 主 河道 时 浅 层 含水 层 的 水 位 国 值 (mm) 8.5 
GW. REVAP.gw 地 下 水 再 蒸发 系数 0.03 
ESCO.hru 土壤 蒸发 补偿 系数 1.008 
ALPHA_BNK.rte 基 流 退 水 常数 (days) 0.45 
SFTMP 降雪 日 的 平均 空气 温度 (Y ) 4.8 
PLAPS 降水 直 减 率 (mm/km) 44.5 
SMFMN 12 月 21 日 的 融雪 因子 (mm/d-%C) 9.407 
SMFMX 6 H 21 日 的 融雪 因子 (mm/d-%C) 0.1302 
TLAPS 温度 直 减 率 ( /km) -4.3039 


CN2:SCS 径流 曲线 值 , Initial SCS runoff curve number for moisture condition I;ALPHA_BF: 基 流 a 因子 ,Baseflow alpha factor;GW_DELAY: 地 下 
水 补给 延迟 时 间 ,Groundwater delay time; GWQMN :地 下 水 汇 入 证 河道 时 浅 层 含水 层 的 水 位 阔 值 , Threshold depth of water in the shallow aquifer 
required for return flow to occur; GW_REVAP ; 4E F ZK-HEz& A AR, Groundwater " revap" coefficient; ESCO: 土壤 蒸发 补偿 系数 , Soil evaporation 
compensation factor; ALPHA, BNK ; 3£ ji; 3H 7K 2$ Zt; Baseflow alpha factor for bank storage; SFTMP ; 降雪 日 的 平均 空气 温度 , Snowfall temperature; 
PLAPS :降水 直 减 率 ,Precipitation lapse rate; $óMFMN; 12 月 21 日 的 融雪 因子 ,Melt factor for snow on December 21; SMFMX:6 月 21 日 的 融雪 因子 ， 
Melt factor for snow on June 21; TLAPS :温度 直 减 率 ,Temperature lapse rate 


3.3 ”模型 评估 

本 研究 选取 (NSE Nash-Sutcliffe Efficiency) 效率 系数 及 R 决定 性 系数 作为 模型 评估 指标 ,以 上 指标 被 广 
泛 的 应 用 于 评估 模型 的 性 能 站 。 其 中 ,Nash-Sutcliffe 效率 系数 是 一 个 正 态 统计 方程 , 它 反 映 了 观测 值 与 对 
应 模拟 值 的 拟 舍 程度 @ NSE 可 利用 方程 (1) HEREDI. 


2.50809 
NSE = 1 -一 — (1) 
2. (Qni a Qna) 
WEF ,万 为 径流 值 变量 , 0, M 0. 分 别 代表 观测 值 及 模式 值 。0, 表示 径流 平均 观测 值 。NSE 变化 范围 为 -oo 至 
1,25 NSE 计算 结果 为 1 时 ,认为 观测 值 与 模拟 值 完全 吻合 ; 当 其 值 为 0.5 一 1 之 间 时 ,模式 结果 可 接受 , 当 NSE 
小 于 0 时 ,认为 模式 结果 较 差 。 
决定 性 系数 RR 表示 变量 之 间 相 关 程 度 ( 见 式 2) : 


[3(09,,-02(0,; -0)] 
ccc 一 - (2) 
2 (070. "X Oeo 


i 
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式 中 , Q, , Q. 和 式 (6) 中 的 含义 相同 ,i 变量 表示 第 i 次 的 模拟 或 实测 值 。 

众多 研究 将 R 20.5 H. NSE>0.5 作为 SWAT 模型 模拟 满意 的 准则 ” , 亦 有 研究 认为 NSE>0.4 也 可 作为 
模型 指标 满意 的 标准 ,本 研究 采用 Moriasi 等 评价 标准 [9 。 即 模型 在 率 定 期 的 月 尺度 模拟 结果 NSE 20.65 或 
日 尺度 模拟 结果 NSE 20.5 ,说明 模型 模拟 结果 可 接受 。 


4 结果 分 析 


41 基于 CMADS+SWAT 模式 的 征 流量 模拟 验证 及 分 析 。 率 定期 
本 研究 利用 CMADS+SWAT 模式 分 别 输出 流域 两 
内 两 个 水 文 控制 站 ( 精 河山 口 控制 站 及 温泉 站 ) 逐 月 径 
流量 ,利用 观测 结果 进行 参数 率 定 与 校准 工作 (图 4 一 
图 7)。 
分 析 发 现 :在 月 尺度 上 ,由 CMADS 驱动 下 的 SWAT 
模型 模拟 结果 在 精 博 河流 域 两 个 控制 站 均 到 达 了 满意 


验证 期 A 95PPU 
A 观测 值 
A 模拟 值 


EN 


o2 


径流 量 Flow out/(m?/s) 


10 20 30 40 50 60 
的 效果 ( 见 表 5)。 由 CMADS 驱动 下 的 模拟 结果 在 精 河 时 间 Time/months 
山口 站 的 NSE 和 R? 均 取得 了 令 人 满意 的 表现 (NSE = 4 CMADS 数据 集 驱动 SA 了 模型 的 精 河山 口 控制 站 逐 月 模 
0.939, R^ 20.942) , MER (2009—2013), dB ,95PPU 代表 预测 不 确定 性 内 监测 数 


在 模型 验证 期 ,虽然 NSE 效率 系数 与 确定 性 系数 EUER 

REER ss HH p ME 也 得 到 了 可 接受 的 效果 ( NSE - 0. Fig.4 The ae and MEN nn d NN 
by SWAT. model driving by CMADS in Jinghe control station, 

904 ,R" 3 0.934) S 相对 精 河 山 H 控制 站 来 次 2 LES 站 在 where the percentage of measured data bracketed by the 95% 

率 定 期 和 验证 期 的 模拟 精度 均 较 前 者 KR, 笔者 认为 , prediction uncertainty ( 95PPU) 

ie Ef eg EL KAKI DE e P AZ ARS] ELE, , EH 

T SWAT 模型 度 日 因子 仅 对 融雪 进行 了 细致 化 考虑 ,导致 流域 在 温 录 站 的 模拟 精度 要 低 于 精 河 控制 站 (冰川 

补给 率 较 温 果 控制 站 小 ) 所 在 子 流域 。 


-" 率 定期 (月 ) TM 验证 期 (月 ) 
a b 

45 y= 0.9976x - 0.0007 40 + y= 1.094x - 0.3648 

40 R? = 0.9423 35 上 R? = 0.9342 
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径流 模拟 值 
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5 精 河山 口 控制 站 月 观测 与 模式 模拟 结果 确定 性 系数 


Fig.5 The certainty factor of the relationship between monthly observation and model simulation data in Jinghe control staion 


当 完 成 月 尺度 率 定 和 验证 后 ,本 研究 将 月 尺度 最 佳 参数 值 带 入 SWAT 模型 进行 日 微调 校准 后 发 现 . 
CMADS+SWAT 模式 径流 模拟 结果 在 两 个 水 文 控制 子 流域 与 日 观测 结果 呈现 很 好 的 一 致 性 。( 兄 网 8 一 图 
11、 表 5)。 

在 模型 率 定期 ,由 CMADS 驱动 SWAT 模式 在 精 河 站 的 日 拟 合 结 果 ( NSE 20.801, R^ 2 0.815) 5i if 5 vi 0 
合 结果 (NSE= 0.796, 尺 =0.791) 也 均 较 好 。 在 验证 期 ,虽然 SWAT 模式 在 两 个 站 点 模拟 结果 均 可 接受 ,而 模 
式 在 精 河山 口 站 模拟 效果 (NSE=0.851 , R^ 20.796) $i 5 Ta tll yi ( NSE- 0.526, R^ 20.592) 表现 更 佳 。 
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4.0 基于 CMADS+SWAT 模式 的 土壤 湿度 及 相关 变量 分 


将 SWAT 模型 进行 精 博 河流 域 本 地 化 后 ,为 从 时 空 
角度 分 析 CMADS+SWAT 模式 在 土壤 湿度 及 融雪 变量 
时 空 演 变 模 拟 能 力 , 同 时 定量 分 析 各 分 量 之 间 响 应 关 
系 , 本 研究 还 以 精 河 子 流域 为 主要 分 析 对 象 ,对 多 个 要 
素 啊 应 关系 进行 了 提取 。 

4.2.4. 精 博 河流 域 融 雪 过 程 与 土壤 湿度 啊 应 分 析 

为 俩 究 精 博 河流 域 融 雪 过 程 对 土壤 湿度 的 影响 ,本 
部 分 提取 精 博 河 整 个 流域 在 2012 年 4 月 8 日 的 土壤 湿 
度 及 其 对 应 融雪 量 空间 变化 图 (图 12) 。 

图 12A 为 当日 土壤 温度 分 布 情况 ,图 12B 为 对 应 
时 期 融雪 空间 分 布 情况 。 为 定量 研究 流域 内 各 地 表 分 
量 间 关系 ,研究 提取 了 精 河 控制 站 子 流域 内 各 地 表 分 量 
(包括 土壤 温度、 潜在 苹 散 发 .降水 量 和 融雪 量 ) 在 2012 
年 4 月 8 日 的 时 间 序 列 分 布 (图 13)。 
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6 CMADS 数据 集 驱 动 SWAT 模型 的 温泉 控制 站 逐 月 模拟 结 
果 (2009 一 2013)， 其 中 ,95PPU 代表 预测 不 确定 性 内 监测 数据 占 
据 的 百分比 

Fig.6 The monthly simulation results simulated by SWAT model 
driving by CMADS in Wenquan control station, where the 
percentage of measured data bracketed by the 95% prediction 
uncertainty ( 95PPU) 


TS n 验证 期 (月 ) 


y= 1.1366x - 0.1751 


q 
wE 
ES 
BP 16 1 
w 3 2 
go 14 
Er , 
9 
e g 
6 4 
6 11 16 21 26 
径流 模拟 值 
Runoff simulated/(m?/s) 
图 7 “温泉 控制 站 月 站 点 观测 与 模式 模拟 结果 确定 性 系数 


Fig.7 The certainty factor of the relationship between monthly observation and model simulation data in Wenquan control station 


通过 对 2012 4E 44H 8 日 精 情 河流 域内 土壤 湿度 分 
布 分 析 发 现 : 流域 整体 填 壤 温度 在 该 日 处 于 侦 湿 状态 
(图 14) 。 其 中 , 温 果 站 控制 子 流域 土壤 湿度 为 18mm 
至 20mm 左右 , 而 精 河 控制 站 及 其 附近 子 流域 土壤 湿度 
达 近 57mm, 遂 过 分 析 流 域内 对 应 融雪 空间 分 布 情况 
(图 14) 后 友 现 ,2012 年 4 月 8 日 , 精 博 河流 域内 出 现 整 
体 较 头 融 必 现 绷 , 而 发 生 融 雪 的 区 域 恰 好 为 西天 山北 坡 
地 区 ,该 区 域 年 内 3 一 4 月 雪 盖 覆盖 丰厚 。 为 定量 分 析 
精 博 河 流域 内 融雪 及 土壤 湿度 量 级 及 其 直接 响应 关系 ， 
本 人 研 究 对 精 河 水 文 站 控制 子 流域 中 土壤 湿度 与 融雪 量 
等 分 量 进 行 了 (以 精 河 山口 控制 站 为 例 ) 时序 量 级 分 析 
人 研究 。 

图 13 给 出 了 精 河 控制 子 流域 内 各 地 表 或 近 地 表 分 
量 ( 洪 在 蒸 散发 ,土壤 湿度 ,降水 量 及 融雪 量 等 ) 在 精 河 
山口 控制 子 流 域 时 序 分 布 状 态 。 分 析 发 现 : 精 博 河流 域 
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8  CMADS 数据 集 驱动 SWAT 模型 的 精 河山 口 控制 站 逐日 模 
拟 结果 (2009 一 2013) ,其 中 ,9SPPU 代表 预测 不 确定 性 内 监测 数 
据 占据 的 百分比 

Fig.8 The daily simulation results ( 2009—2013) simulated by 
SWAT model driving by CMADS in Jinghe control station, where 
the percentage of measured data bracketed by the 95% prediction 
uncertainty ( 95PPU) 
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21 期 孟 现 勇 ” 等 :基于 CMDAS 驱动 SWAT 模式 的 
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Fig.9 The certainty factor of the relationship between daily observation and model simulation data in Jinghe.control staion 


H 2012 年 3 月 中 旬 开 始 出 现 融雪 现象 ,并 于 2012 年 4 
月 8 日 出 现 较 大 融雪 现象 。 其 中 , 精 河 站 控制 子 流域 融 
雪 量 达 18mm/d。 当 日 精 河 控制 站 子 流域 土壤 湿度 也 
达到 其 年 内 较 高 水 平 ( 约 56mm) 。 分 析 表 明 : 融 雪 期 
内 , 精 博 河流 域 土壤 湿度 增加 的 最 主要 贡献 来 自 于 融 
雪 , 仅 极 少 量 来 自 降 水 。 此 外 ,土壤 湿度 与 潜在 蒸发 数 
值 上 呈 反 比 关系 。 
4.2.2，” 精 博 河流 域 降水 与 土壤 湿度 响应 分 析 

前 文 提 到 精 河 站 控制 子 流 域 融雪 过 程 对 流域 内 填 
培 湿 度 具 有 较 大 影响 。 笔 者 认 为 降水 量 对 土壤 湿度 页 
献 亦 不 容 忽视 。 笔 者 在 此 重点 分 析 了 夏季 与 秋 来 降水 
对 土壤 湿度 的 影响 。 由 于 精 博 河流 域 土 壤 湿 度 在 融雪 
后 期 至 年 内 9 月 会 出 现 较 长 时 期 波动 ,同时 10 明 中 后 
期 将 发 生 陡 然 上 升 现象 并 保持 恒定 等 现象 ) 笔 者 针对 以 
上 两 阶段 降水 量 与 土壤 湿度 响应 关系 进行 时 空 提 取 分 
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10 CMADS 数据 集 驱 动 SWAT 模型 的 温泉 控制 站 逐日 模拟 
结果 (2009 一 2013) ， 其 中 ,9SPPU 代表 预测 不 确定 性 内 监测 数据 
占据 的 百分比 

Fig.10 The daily simulation results ( 2009 一 2013 ) simulated by 
SWAT model driving by CMADS in Wenquan control station, 
Where the percentage of measured data bracketed by the 95% 
prediction uncertainty ( 95PPU) 
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Fig.11 The certainty factor of the relationship between daily observation and model simulation data in Wenquan control station 


析 ( 图 14). 2g Scu E LU LET T5 ri] ns E PRESE 2E LEER P RARE Tt CH : UR dC SS DR) 演变 情况 ,研究 还 


提取 了 精 河 山口 控制 站 冻 土 和 雪 深 逐日 台 站 观测 资料 ( 
图 14 给 出 了 精 博 河流 域 土壤 湿度 与 降水 量 在 日 内 
况 。 其 中 ,图 14 为 2010 年 6 月 22 日 精 博 河 流域 土壤 


2010 至 2011 4E) (图 16) 作 为 验证 分 析 。 
(2010 年 6 月 22 日 和 2011 年 10 月 21 日 ) 空 间 分 布 情 
湿度 分 布 情况 ,图 14 为 对 应 流域 内 降水 量 分 布 情况 。 


分 析 发 现 :当日 精 河山 口 控制 站 子 流 域 土壤 湿度 达 24.8mm, 精 河山 口 控制 站 子 流 域 的 降水 量 达 30.1mm。 分 
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12 ” 精 博 河 流域 土壤 湿度 与 融雪 量 空间 相关 分 析 


Fig.12 The correlation analysis of soil moisture and snow melting in Jinghe River Basin 
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13， 精 博 河流 域 各 类 分 量 时 间 序列 分 析 ( 以 精 河山 口 控制 站 为 例 ) 


Fig.13 The correlation analysis in Jinghe River Basin (taking the control station of Jinghe as an example) 
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14 精 博 河流 域 土壤 湿度 与 降水 量 空间 相关 分 析 


Fig.14 The correlation analysis of soil moisture and precipitation in Jinghe River Basin 


析 表 明 :降水 量 在 日 内 大 幅度 增加 导致 该 子 流域 土壤 湿度 在 短 时 期 快速 攀升 。 类 似 情 况 同 时 出 现 于 2011 4E, 
图 14 给 出 了 精 博 河 流域 2011 年 10 月 21 日 土壤 湿度 的 空间 分 布 状态 ,图 14 为 对 应 时 期 降水 量 在 整个 流域 
的 分 布 状态 。 该 时 有 段 精 河 控 制 站 及 温 果 控制 站 所 在 的 子 流域 土壤 湿度 较 大 ,其 中 温 果 站 控制 站 所 在 子 流域 土 
HERBES 18mm 左右 , 精 河山 口 控 制 站 所 在 子 流域 土壤 湿度 达到 近 22mm。 分 析 还 发 现 , 在 流域 内 土壤 湿度 达 
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到 该 时 期 较 高 水 平 的 同时 ,降水 过 程 在 精 博 河 整 个 子 流 域 均 有 发 生 。 为 准确 的 验证 土壤 湿度 与 降水 量 在 量 级 
上 的 响应 关系 ,本 人 研究 提取 了 精 河 山口 控制 站 (2010 年 至 2012 年 ) 其 他 地 表 分 量 ( 湾 在 蒸 散发 土壤 湿度 LE 
水 量 及 融雪 量 ) 对 应 时 序 关系 曲线 (图 15) 。 


表 5 CMDAS+SWAT 模式 径流 输出 的 月 .日 拟 合 结果 评估 
Table 5 The assessment under monthly and daily scale simulated by CMADS+SWAT mode 


时 间 尺 度 精 河 山口 控制 站 Jinghe control station 温 果 控制 站 Wenquan control station 

Time scale 率 定期 NS 率 定 期 验证 期 NS 验证 期 率 定 期 NS 率 定 期 验证 期 NS 验证 期 
Calibration NSE Calibration R? Validation NSE Validation R? Calibration NSE Calibration R? Validation NSE Validation R? 

H Month 0.939 0.942 0.904 0.934 0.917 0.914 0.659 0.844 

H Day 0.801 0.815 0.802 0.851 0.796 0.791 0.526 0:592 


NSE : 纳什 效率 系数 ，Nash-Sutcliffe Efficiency 


分 析 发 现 ,在 2010 年 6 月 22 日 ,流域 内 精 河 山口 控制 站 子 流 域 降 水 量 接近 30mm , 流域 内 土壤 湿度 亦 擎 
JF 28 24.8mm( KI 15) , 2011 年 10 H 21 日, 精 河山 口 控制 站 所 在 子 流域 内 降水 量 达 近 和 mm ,于 此 同时 , 子 流 
域内 土壤 湿度 值 也 快速 攀升 至 22mm 左右 ,土壤 湿度 在 此 后 保持 恒定 状态 (图 15) 。j 秋 未 时 期 , 精 河 山口 控制 
子 流域 土壤 湿度 快速 上 升 现 象 充分 说 明 ,年 末 降 水 对 土壤 湿度 变化 贡献 较 类 。 
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15 ” 精 博 河流 域 土壤 湿度 与 降水 量 时间 序 列 分 析 ( 以 精 河 山口 控制 站 为 例 ) 
Fig.15 The correlation analysis of Soil moisture and precipitation in Jinghe River Basin (taking the control station of Jinghe as an 


example) 


以 年 内 循环 角度 对 土壤 湿度 变化 进行 分 析 后 发 现 :人 研究 区 内 每 年 3 至 4 月 为 融雪 期 ,流域 内 雪 深 经 历 了 
先 上 升 后 下 降 的 过 程 , 并 在 融雪 期 结束 时 雪 深 接近 0m, 这 与 本 研究 在 精 河 站 提取 的 雪 深 观测 数据 一 致 (图 
16) 。 融 雪 期 后 期 气候 转 暖 ,此 时 多 十 期 来 临 ,造成 研究 区 内 土壤 泥 度 快速 葡 升 ,由 于 空气 温度 在 该 时 段 较 
高 , 冻 土 发 生 融化 现象 (图 16) 。 共 发 同时 增加 并 加 速 土壤 水 减少 ,导致 该 时 期 土壤 湿度 呈现 波动 变化 趋势 。 
10 月 后 期 ,流域 由 于 冷 湿 气流 过 境 出 现 大 范围 降水 (降雪 ) ,此 时 土壤 含水 量 将 大 幅 增 加 ,伴随 冷 空 气 来 临 而 
促使 蒸发 量 大 幅度 下 降 。 此 时 ,土壤 液态 含水 量 冻 结 ( 每 年 11 月 至 12 月 ) 并 产生 冻 土 。 季 闻 性 冻 土 将 在 来 
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16 2010 年 .2011 年 精 河 控制 站 雪 深 、 冻 土 观测 值 
Fig.16 The snow depth and frozen soil depth observation in Jinghe control station ( Year 2010-and 2011) 


年 1 一 2 月 到 达 峰 值 直 至 次 年 融雪 期 再 次 陡然 下 降 , 冻 土 再 次 融化 ,A 此 时 ,一 部 份 冻 主将 转换 变 成 土壤 液态 含 
水 量 。 

综 上 所 述 ,流域 内 土壤 湿度 在 每 年 3 一 4 月 份 到 达 其 第 一 次 较 电 水 平 ,这 主要 是 由 于 流域 内 积 雪 、 冻 土 融 
化 而 造成 的 。 融 雪 期 结束 后 ,由 于 降水 量 增多 ,空气 温度 上 天 导致 主 壤 湿度 呈现 上 下 波动 态势 ,直至 10 月 中 
旬 冷 空气 过 境 产 生 较 大 降水 ( 雪 ) 并 最 终 使 土壤 水 转换 为 还 塘 , 至 次 年 融雪 期 ,土壤 液态 水 再 次 增加 直至 融雪 
期 结束 。 


s ”讨论 与 结论 


本 文 利用 CMADS 驱动 SWAT 模型 对 精 精 河 法 域内 两 个 水 文 控制 站 进行 了 径流 精度 模拟 验证 工作 ,利用 
本 地 化 的 SWAT 模型 分 析 了 精 河 控制 站 所 在 子 流域 各 类 地 表 分 量 (包括 土壤 湿度 \ 洪 在 区 散发 .降水 量 和 融 
雪 量 ) 时 空 变化 规律 ,加 深 了 对 该 巴 流 域内 名 分 量 要 系 的 空间 演变 规律 ,解释 了 人 研究 区 内 土壤 湿度 在 不 同 季 
节 的 变化 机 制 。 本 人 研究 后 期 将 继续 考虑 将 更 多 观测 数据 (如 卫星 遥感 数据 ) 同 化 到 SW AT 模型 输出 数据 中 ， 
以 提高 SWAT 模式 的 结果 可 年 性 蛙 在 现 有 阶段 ,利用 CMADS 数据 集 驱 动 SW AT 模型 对 研究 高 分 辨 率 气 象 数 
据 驱 动 复杂 水 文 模型 的 研究 意义 较 大 ,预计 该 数据 集 将 为 我 国 的 水 文生 态 细致 化 模拟 研究 提供 数据 文 掉 。 


参考 文献 ( References) : 


[1] WAR, mhi, EH f, 刘 志 辉 , 吉 晓 楠 , 余 丹 林 . Xr 60 年 气候 变化 及 人 类 活动 对 艾 比 湖 流 域 水 资源 的 影响 . 水 文 , 2015, 35(2) : 
90-96. 

[2] Bere A A, Famiglietti J S, Walker J P, Houser P R. Impact of bias correction to reanalysis products on simulations of North American soil moisture 
and hydrological fluxes. Journal of Geophysical Research, 2003, 108( D16) : 4490. 

[3] Maurer E P, Wood A W, Adam J C, Lettenmaier D P, Nijssen B. A long-term hydrologically based dataset of land Surface fluxes and states for the 
conterminous United States. Journal of Climate, 2002, 15(22) : 3237-3251. 

[4] Fekete B M, Vórósmaety C J, Roads J O, Willmott CJ. Uncertainties in precipitation and their impacts on runoff estimates. Journal of Climate, 
2004, 17(2) : 294-304. 

[5] Sheffield J, Ziegler A D, Wood E F, Chen Y B. Correction of the high-latitude rain day anomaly in the NCEP-NCAR reanalysis for land surface 
hydrological modeling. Journal of Climate, 2004, 17(19) : 3814-3828. 

[6] Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L, Iredell M,Saha S, White G, Woollen J, Zhu Y, Leetmaa A, Reynolds B, 
Chelliah M, Ebisuzaki W, Higgins W, Janowiak J, Mo K C, Ropelewski C, Wang J, Jenne R, Joseph D. The NCEP/NCAR 40-year reanalysis 
project. Bulletin of the American Meteorological Society, 1996, 77(3) : 437-472. 


http ://www.ecologica.cn 


2] 期 孟 现 勇 ” 等 :基于 CMDAS 驱动 SWAT 模式 的 精 博 河 流域 水 文 相 关 分 量 模拟 验证 及 分 析 7127 


We 


[8] 


RA 


[10] 


[11] 


[12] 
[13] 
[14] 
[15] 


[16] 


[17] 


[18] 


[19] 
[20] 


[21] 


[22] 
[23] 
[24] 
[25] 


[26] 


[27] 


[28] 
(29 | 


[30] 


[21] 


Kanamitsu M, Ebisuzaki W, Woollen J, Yang S K, Hnilo J J, Fiorino M, Potter G L. NCEP-DEO AMIP-II reanalysis ( R- 2) . Bulletin of the 
American Meteorological Society, 2002, 83(11) ; 1631-1643. 

Uppala S M, KÁllberg P W, Simmons A J, Andrae U, Da Costa Bechtold V, Fiorino M, Gibson J K, Haseler J, Hernandez A, Kelly G A, Li X, 
Onogi K, Saarinen S, Sokka N, Allan R P, Andersson E, Arpe K, Balmaseda M A, BeljaarsA C M, Van De Berg L, Bidlot J, Bormann N, 
Caires S, Chevallier F, Dethof A, Dragosavac M, Fisher M, Fuentes M, Hagemann S, Hólm E, Hoskins B J, Isaksen L, Janssen P A E M, Jenne 
R, Mcnally A P, Mahfouf J F, Morcrette J J, Rayner N A, Saunders R W, Simon P, Sterl A, Trenberth K E, Untch A, Vasiljevic D, Viterbo P, 
Woollen J. The ERA-40 re-analysis. Quarterly Journal of the Royal Meteorological Society, 2005, 131 (612) : 2961-3012. 

Onogi K, Tsutsui J, Koide H, Sakamoto M, Kobayashi S, Hatsushika H, Matsumoto T, Yamazaki N, Kamahori H, Takahashi K, Kadokura S, 
Wada K, Kato K, Oyama R, Ose T, Mannoji N, Taira R. The JRA-25 reanalysis. Journal of the Meteorological Society of Japan, 2007, 85(3) : 
369-432. 

TR, 师 春香 , 刘 时 银 , 王 浩 , 雷 晓 辉 , 刘 志 辉 , 吉 晓 楠 , 蔡 思 宇 , 赵 求 东 . CMADS 数据 集 及 其 在 流域 水 文 模型 中 的 驱动 作用 一 使 以 黑河 
流域 为 例 . 人 民 珠 江 , 2016, 37(7) : 1-19. 

Pal J S, Eltahir E A. Pathways relating soil moisture conditions to future summer rainfall within a model of the land-atmosphere system. Journal of 
Climate, 2001,14(6) ; 1227-1242. 

赵 天 保 , 艾 丽 坤 , 冯 钊 明 . NCEP 再 分 析 资 料 和 中 国 站 点 观测 资料 的 分 析 与 比较 . 气候 与 环境 研究 42004, 9(2):« 278-294. 

赵 天 保 , FRA. 几 种 再 分 析 地 表 气 温 资料 在 中 国 区 域 的 适用 性 评估 . 高 原 气 象 , 2009, 28(3): 595-607. 

赵 天 保 , 符 深 斌 . 中 国 区 域 ERA-40, NCEP-2 再 分 析 资 料 与 观测 资料 的 初步 比较 与 分 析 . 氢 候 与 环境 人 研究，2006，11(1) : 14-32. 

施 晓 晖 , REEE, 谢 立 安 . NCEP/NCAR 再 分 析 风 速 .表面 气温 忠平 在 中 国 区 域 气候 变化 研究 中 的 如 信和 度 分 析 . 气象 学 报 , 2006, 64(6): 
709-722. 


Fischer G, Nachtergaele F, Prieler S, van Velthuizen H T, Verelst L, WibergD:2008. Global agro-ecological zones assessment for agriculture 
( GAEZ 2008). Laxenburg, Austria, Rome, Italy; MASA, FAO. 

刘 时 银 , WIRE, 郭 万 钦 , VERON], LEZHE, BRE, MBE, RUR. 基于 第 二 次 冰川 编目 的 中 国 冰川 现状 . 地 理学 报 , 2015, 70 
(1) : 3-16. 

Monteith J L. 1965. Evaporation and environment. The state and movement of water in living organisms//The 19th Symposium of the Society on 
Experimental Biology. Swansea; Cambridge University Press. 

Allen R G. A Penman for all seasons. Journal of Irrigation.and Drainage Engineering, 1986, 112(4) : 348-368. 

Allen R G, Jensen ME, Wright JL, Burman R D. Operational estimates of reference evapotranspiration. Agronomy Journal, 1989, 81(5): 
650-662. 

Priestley C H B, Taylor R J. On the assessment 6f surface heat flux and evaporation using large-scale parameters. Monthly Weather Review, 1972, 
100(2) : 81-92. 

Hargreaves G H. Moisture availability and crop production. Transactions of the ASAE, 1975, 18(5) : 980-984. 

Hargreaves G H, Samani Z A: Reference crop evapotranspiration from temperature. Applied Engineering in Agriculture, 1985, 1(2) : 96-99. 
Hargreaves G H, Samani Z A. Estimating potential evapotranspiration. Journal of the Irrigation and Drainage Division, 1982, 108(3) : 225-230. 
Hargreaves G L, Hargreaves GH, Riley J P. Agricultural benefits for Senegal River Basin. Journal of Irrigation and Drainage Engineering, 1985, 
111(2) : 113-124. 

Abbaspour-K;- Yang J, Reichert P, Vejdani M, Haghighat S, Srinivasan R. SWAT-CUP. SWAT calibration and uncertainty programs, Swiss 
Federal Institute of Aquatic Science and Technology ( EAWAG) , 2008, Zurich, Switzerland. 

Nash J E, Suteliffe R M. River flow forecasting through conceptual models part 1-a discussion of principles. Journal of Hydrology, 1970, 10(3) : 
282-290. 

Schaefli B, Gupta H V. Do Nash values have value? Hydrological Processes, 2007, 21 (15): 2075-2080. 

Santhi C, Arnold J G, Williams J R, Hauck L M, Dugas W A. Application of a watershed model to evaluate management effects on point and 
nonpoint source pollution. Transactions of the ASAE, 2011, 44(6) : 1559-1570. 

Ahmad H M N, Sinclair A, Jamieson R, Madani A, Hebb D, Havard P, Yiridoe E K. Modeling sediment and nitrogen export from a rural 
watershed in Eastern Canada using the soil and water assessment tool. Journal of Environmental Quality, 2011, 40(4) : 1182-1194. 

Moriasi D N, Arnold J G, Van Liew M W, Bingner R L, Harmel R D, Veith T L. Model evaluation guidelines for systematic quantification of 
accuracy in watershed simulations. Transactions of the ASABE, 2007, 50(3): 885-900. 


http ://www.ecologica.cn 


